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SHOCK-WAVE SHAPES AND SURFACE PRESSURES ON HEMISPHERE- 


AND CONE-CYLINDERS IN LOW-DENSITY HYPERSONIC 
FLOW INCLUDING ABLATION 
Daryl J. Monson and Donald M. Kuehn 
Ames Research Center 


SUMMARY 


Shock-wave shapes and surface pressures were determined in air for highly 
cooled hemisphere- and cone-cylinders with and without ablating noses. The 
test Mach number was 14. The test Reynolds numbers of 9000 and 1000, based on 
the cylinder diameter, were estimated to be in the boundary-layer and merged- 
layer flow regimes, respectively. Comparisons with the predictions of 
viscous-interaction theory have been made and the range of validity for the 
theory has been assessed. 

The results show that in the merged-layer regime viscous -interaction 
theory can predict the shock-wave shapes for these bodies but fails to pre- 
dict the surface pressures on the bodies. With merging, the surface pressures 
on the hemisphere-cylinder decreased below the predicted inviscid pressures. 
The behavior of the shock waves and surface pressures for these bodies in the 
merged-layer regime is qualitatively the same as that previously observed for 
sharp bodies . 

Nose ablation displaced the shock waves outward and increased the surface 
pressures on the bodies for a short distance behind the noses. These trends 
were qualitatively predicted by viscous-interaction theory. 


INTRODUCTION 


The flow over sharp flat plates, wedges, and cones successively enters 
the boundary- layer, weak- and strong-interaction, and merged-layer regimes as 
the rarefaction of the flow increases. Viscous -interaction theories (e.g., 
ref. 1) have been shown to predict the shock-wave locations on these bodies 
in all the regimes (ref. 2) but to overpredict the measured surface pressures 
in the merged-layer regime (e.g., refs. 2, 3, and 4). McCroskey, et al . 

(ref. 5) , attribute this failure to a departure from a Rankine-Hugoniot 
shock wave when the boundary layer merges with the shock wave. 

Probstein (ref. 6) classified the flow over blunt bodies into two 
regimes. The first includes the boundary-layer, vorticity-interaction, and 
viscous-layer regimes; the second includes the incipient merged-layer and 
fully merged-layer regimes. In this report they shall be referred to as the 



boundary- layer and merged-layer regimes, respectively. The usual viscous- 
interaction theory for blunt bodies numerically calculates the boundary- layer 
displacement thickness over a body, adds this to the body thickness, and cal- 
culates the inviscid pressures on the "new" body shape. The theory has suc- 
cessfully predicted surface pressures on a variety of blunt bodies in the 
boundary- layer regime (e.g., refs. 7, 8, and 9) . There are indications that 
the theory may fail at low densities, however. Marchand, et al . (ref. 7) , 
substantially overpredicted pressures on a highly cooled, slender blunt cone 
that may have been in the merged-layer regime, and were unable to offer an 
explanation for the disagreement. Kuehn (ref. 10) measured pressures that 
agreed with inviscid theory for highly cooled hemisphere- and cone-cylinders 
at conditions that may have been merged. He attributed the lack of viscous 
induced pressures to wall cooling which resulted in a negligible displacement 
thickness of the boundary layer over the bodies. Thus, at least for surface 
pressures, viscous-interaction theories appear to fail for blunt bodies at 
low densities just as they do for sharp bodies. Questions remain, however, 
on when and for what reason this failure will occur, and whether it will also 
occur for other properties such as shock-wave shapes. 

In this report, the results of an experimental investigation of the 
shock-wave shapes and surface pressures for ablating and highly cooled non- 
ablating hemisphere- and cone-cylinders in low-density hypersonic flow are 
presented and compared with predictions made by simple viscous-interaction 
theories. The data were obtained as part of a larger study of flow separation 
ahead of flares on a cylinder (ref. 11) . In reference 11, the attached-flow 
pressure data were used in determining flow separation and were not analyzed 
in themselves, as is done in the present study. The experiments were con- 
ducted at a Mach number of 14 and at Reynolds numbers of 9000 and 1000 based 
on the cylinder diameter. The purpose of the report is threefold: to present 
measured surface pressures and shock-wave shapes for two highly cooled blunt 
bodies at conditions in the boundary- layer and merged-layer regimes; to com- 
pare the results with predictions from viscous -interaction theory, and examine 
under what conditions and for what reasons the theory fails; and to give some 
preliminary indications of the effects of nose ablation on the surface 
pressures and shock-wave shapes. 


EXPERIMENTAL TECHNIQUE 


Test Facility and Conditions 

The experiments were performed in the Ames Entry Aerodynamics Tunnel 
sketched in figure 1 . The nominal test Mach number, based on pitot-pressure 
surveys and nonequilibrium nozzle calculations (ref. 12), was 14. The test- 
stream calibration showed a centerline Mach number gradient of 1.8x10 2 per 
cm. Additional facility and stream calibration details are given in refer- 
ence 11. The present investigation was conducted at two conditions: (1) at 

a reservoir pressure of 67 atmospheres, a total enthalpy of 4.2xl0 6 J/kg, and 
a Reynolds number. Re,* d , of 9000 ’ and ^ at a reservoir pressure of 10 
atmospheres, a total enthalpy of 8.2x10® J/kg, and a Reynolds number, 
of 1000. 
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Models 


Shaded areas -ablating surfaces 



radius = 0.033 d d = 3.8 1 cm 

Cone- cylinder 

Figure 2.- Test models. 


The model configurations used in 
these tests are shown in figure 2. The 
cylinder was thick -walled stainless 
steel, the nonablating noses were cop- 
per, and the ablating noses were 
Teflon, Delrin, and high-density pheno- 
lic nylon. The static-pressure ori- 
fices (dia. = 0.24 cm) were located 
along the cylinder with a longitudinal 
spacing of 0.57 cm. Thin-walled, 0.46 
cm i.d. stainless-steel tubing connected 
the orifices to pressure cells located 
2 meters from the model. The models 
were protected by a retractable shield 
while the tunnel was started. 


Surface Pressure, Shock-Wave Shape, and 
Ablation-Rate Measurements 

Surface pressures were measured with capacitance-type pressure cells 
to an accuracy of 5 percent at the lowest pressures measured (30 p Hg) and to 
1 percent at the highest pressures (1000 p Hg) . Pre-run out-gassing procedures 
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eliminated pressure gradients in the lines, and during a run, the pressures 
were allowed to reach equilibrium (typically within 20-40 sec) before data 
were taken. The measured pressures were corrected for longitudinal test- 
stream gradients (ref. 13) and low-density orifice effects (ref. 14). The 
orifice corrections ranged from 0 percent near the shoulder of the 
hemisphere-cylinder to 77 percent far downstream on the cone-cylinder at 
Re^ d = 1000, and resulted in scaling the measured pressures upward. The 
latter corrections were made using wall heat-transfer rates calculated by a 
boundary- layer program to be discussed later. For the present highly cooled 
bodies, the corrections are insensitive to the heat transfer. The same 
boundary- layer calculations predicted surface shear stresses that were, in 
all cases, less than the local pressure, thus allowing Potter's method to be 
applied. The low-density effects on stagnation pressure (ref. 15) and the 
effects of temperature gradients along the pressure lines (ref. 16) were 
negligible. 

Shock-wave shapes were measured from photographs of the natural 
luminescent flow field around the models during the runs. The photographs 
are not included in this report because the shock waves were not visible in 
the reproductions. This method was used because normal visualization techni- 
ques, such as schlieren, are not adequate at the low densities encountered. 
Pitot-pressure surveys verified that the shock wave over the cylindrical por- 
tion of the models was located where the luminescence changed intensity. The 
measured shock -wave shapes were corrected for the longitudinal test-stream 
gradients (ref. 13) . 

Ablation rates for the ablating noses were obtained by weighing the 
models before and after a run and dividing the difference by the run time 
corrected for the estimated starting time for ablation. Quoted rates are 
averages of these values, usually over several runs of various durations. 

The rates should be considered as only approximate since the materials may 
not have ablated at a constant rate during the runs, and the phenolic nylon 
and Delrin were observed to lose an unknown fraction of their mass by melt. 

To achieve ablation rates as close to steady state as possible, data for the 
ablating bodies were taken after about 30 seconds of run time. In some 
cases, the Delrin noses were blunted considerably by the end of a run. 


RESULTS AND DISCUSSION 


Flow Regimes 

The flow regimes over the noses of the test models may be classified 
according to the criterion of Ahouse and Harbour (ref. 17) that merging for 
blunt bodies occurs for V > 0.3. For the hemisphere-cylinder model, the 
lower Reynolds number test condition corresponds to V = 0.47; the other to 
V = 0.17. This clearly puts the nose of this configuration in the boundary- 
layer regime at Re^d = 9000 and in the merged-layer regime at Re^d = 1000. 
Application of the above criterion to the cone-cylinder model indicates that 
the blunt tip is merged at both conditions. Classification of the flow 
downstream of the noses is discussed below. 
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The just described classifications are consistent with the behavior of 
the ablation vapors that were observed during testing (see photographs in 
ref. 11) . The vapors should diffuse approximately as far out as the boundary- 
layer edge on a body (ref. 18). At Re^ ^ = 9000, the vapors filled one-third 
to one-half of the shock layer over the cylinder for both bodies. At 
Re^ d = 1000, the vapors completely filled the shock layer on both bodies as 
far back as two or three cylinder diameters from the nose-cylinder junction. 
Thus, it is concluded that both bodies are in the boundary- layer regime at 
Re Mj( j = 9000 and in the merged- layer regime at Re^d = 1000. 


\ Shock -Wave Shapes 

The measured shock-wave shape data are presented in figures 3 and 4, 
together with the predicted inviscid and viscous induced shock shapes. Mea- 
surements were not always possible near the nose because the film contrast 
was insufficient to show the shock waves. For some cases, this may have been 
due to merging of the shock wave and viscous layer resulting in a "shock" that 
is thick and diffuse. 

The inviscid shock waves were obtained with a perfect-gas blunt-body and 
characteristics program of Inouye, et al . (ref. 19). The perfect-gas assump- 
tion was justified by nonequilibrium blunt-body calculations derived by a 
program of Garr and Marrone (ref. 20) , which predicted frozen shock layers 
around the bodies at both test conditions . 

The viscous induced shock wave for each model was obtained by calculating 
boundary- layer displacement thickness over the nose using a perfect-gas 
boundary- layer program of Clutter and Smith (ref. 21), adding this thickness 
to the geometric nose thickness, and calculating the inviscid flow field 
around the "new" nose. The boundary layer on the cylinder cannot affect the 
shock-wave shape in the region of interest since characteristic waves origi- 
nating on the cylinder cannot reach the shock wave for several diameters down- 
stream of the nose for either body. The edge conditions used in the boundary- 
layer program were the inviscid wall values. Refinements in these conditions 
were not felt to be justified since the application of boundary- layer theory 
to the present low densities is questionable. To calculate boundary layers 
for the ablating noses the total measured mass-injection rates were matched 
and the distributions of mass injection and surface temperature were estimated 
from the predicted nonablating heating-rate distributions and materials data 
from references 22 and 23, and unpublished data from Winovich. The molecular 
weights of the ablation vapors were not matched, since the program assumed 
air-into-air injection. (From the above references, the molecular weights of 
Teflon, Delrin, and phenolic nylon are 100, 15, and 18, respectively. Since 
the molecular weight of Teflon differs considerably from that for air, calcu- 
lations are not presented for it in this report.) Displacement thicknesses 
for the ablating cases were calculated according to the definition of Hayasi 
(ref. 24) for boundary layers with blowing. 


The nonablating shock-wave shape data in figure 3 show that the shock 
wave for the hemisphere-cylinder is independent of Reynolds number for both 
test conditions, whereas the shock wave for the cone-cylinder is displaced 
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outward at both conditions with the greater amount occurring at the lowest 
Reynolds number. These effects are consistent with the boundary- layer calcu- 
lations which show that the displacement thickness over the cooled hemisphere 
is always negligible, but on the cone is significant at both conditions and 
largest at the lower Reynolds number. The calculated shock -wave shapes agree 
well with the measured results. The agreement with theory into the merged- 
layer regime is similar to that previously observed for sharp flat plates and 
wedges . 


Measured 

Measured 

O Reao.d = 9000 

O Reco d = 9000 

D Reco,d = 1000 

□ Re® j - 1 000 

Calculated 

Calculated 

[inviscid 

Inviscid 


Reco ( d ~ 9000 

Reco <1 = 1 000 

Reco.d = ,000 




Figure 3.- The shock shapes for the nonablating bodies; = 14. 


The shock-wave-shape data in figure 4 show that for ablating noses that 
ablation causes an outward movement of the shock wave. This is understandable 
since calculations show that ablation increases the boundary-layer displace- 
ment thickness on the nose for all cases. The extent of displacement varies 
considerably, however, and may be influenced by Reynolds number, ablation 
rate, ablation-gas molecular weight, and surface temperature. Since these all 
vary together, their individual effects cannot be separated from the data. 

The viscous-interaction calculations, which include the effects of all the 
variables except molecular weight, are observed to predict the shock positions 
adequately for the phenolic-nylon noses at both test conditions. Calculations 
were not made for the other ablation materials. 


Surface Pressures 

The measured surface-pressure data are presented in figures 5 and 6 
along with the predicted inviscid and viscous induced pressures. Measurements 
were not always obtained near the shoulder since they were not needed in the 
original separation study in which the present data were obtained. Inviscid 
pressures and viscous induced pressures labelled "present method" were 
obtained as previously discussed for the theoretical shock -wave shapes by 
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Measured 
O Nonablating 

□ Teflon ablating, — — — r = O.IO 
P<X) U C0 A 

O Delrin ablating, — — r- = 0.17 

pco u co A 

A Phenolic nylon ablating, — — r- = 0.04 

ptD»(D* 


Calculated 
Inviscid 



Measured 
O Nonablating 

□ Teflon ablating, — — -r = 0.01 

Pco Uco A 

O Delrin ablating, — — t - = 0.08 

Pco^cd* 

A Phenolic nylon ablating, — — t- = 0.04 

Pco u co A 


Calculated 


Inviscid 

m 

PcD U C0 


A 


= 0.04 



(a) Hemisphere-cylinder, Re m d = 9000. 


(c) Cone-cylinder, Re m d = 9000. 


Measured 
O Nonablating 

□ Teflon ablating, — — — t = 0.24 
Pco u co A 

O Delrin ablating, — — — * = 0.44 

pco uco A 

A Phenolic nylon ablating, — — — t- = 0.18 
Pco u aD A 


Measured 
O Nonablating 

□ Teflon ablating, — 7 } — -r = 0.01 
Pco 

O Delrin ablating, — — t- *0.22 

Pco u co A 

A Phenolic nylon ablating, „ ™ — r =0.11 
Pay u co a 


Calculated 


Calculated 


Inviscid Inviscid 



Figure 4.- The shock shapes for the bodies with ablating noses; = 14. 
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extending the calculations onto the cylinder. This procedure was modified 
for ablating noses, as will be discussed later. This method was not applied 
to the cone-cylinder, since according to pictures of the flow of ablation 
vapors over the corner, using inviscid-wal 1 edge conditions in the boundary- 
layer calculations resulted in what appeared to be unreasonable growth of the 
boundary layer behind the nose-cylinder junction. A second viscous- 
interaction theory was applied to the nonablating bodies in a modified flat- 
plate analysis that Wagner and Watson (ref. 25) used successfully for 
predicting the measured pressures on adiabatic-walled cone-cylinders. 


The nonablating pressure data in figure 5 agree well with the 
vis cous - interaction calculations at R6co 5 d = 9000 (except, perhaps, very near 
the nose), but are significantly below the viscous -interaction predictions at 
Rea, d = 1000. In fact, the pressures on the hemisphere-cylinder at 
Reo/d = 1000 are even lower than the theoretical inviscid pressures. This 
failure of theory to predict surface pressures in the merged-layer regime is 
similar to a trend previously observed for flat plates, wedges, and cones. 

The predicted increase of surface pressures at Re OT} d = 1°00 is due to the 
calculated increased shock-layer displacement caused by the boundary- layer 
growth at that condition. Increased shock- layer displacement, at least over 
the nose of the cone-cylinder, was verified experimentally by the previous 
results which showed increased shock-wave displacement for that body at 
Re^ d = 1000 (see fig. 3(b)). Thus, the failure of viscous -interaction theory 
to°°predict surface pressures at this condition does not appear to be explain- 
able by a decrease in shock-layer displacement on the bodies, as was proposed 
previously by Kuehn (ref. 10) to explain his data. The explanation for the 
failure is logically the same as has been verified for sharp bodies; that is, 
the departure from a Rankine-Hugoniot shock wave upon the merging of the 
shock wave and the boundary layer within the shock layer. This phenomenon 
would also explain why Marchand, et al . (ref. 7), measured pressures on blunt 
cones that were independent of Reynolds number in low-density flow. 

Measured 

O Reco.d = 9000 
□ Reco,d=*000 


Calculated 
Inviscid 

Reco.d = "00- present method 


Fteco^ " 3000, ref. 25 

Recod = '000, present method 

Recod = 1000, ref. 25 



(a) Hemisphere-cylinder. 


Measured 

O R e co,d = 9000 

O Reco d “ 1 000 

Calculated 
Inviscid 


Reco td = 9000, ref. 25 

Reco d " 1000, ref. 25 



(b) Cone-cylinder. 


Figure 5.- The surface pressures for the nonablating bodies; M*, = 14. 
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Surface pressure was not measured near the nose of most of the ablating 
models in figure 6; however, some limited conclusions can be made from avail- 
able measurements. The predicted viscous induced pressures are presented for 
some representative hemisphere -cylinders (no method was available for predict- 
ing pressures on the cone-cylinder) . Since the viscous-interaction theories 
were shown to fail for the nonablating bodies at Re^d = 1000, the theory 
used here is modified to illustrate better the effects of ablation predicted 
when the difference between the predicted ablating and nonablating pressures 
is added to the measured nonablating pressures. The theoretical inviscid 
pressures are also shown on each figure. 

Measured 


O Nonablating 

□ Teflon ablating, — — — t- = O.IO 
Pco u ao A 

O Oelrin ablating, — — t ■ =0.17 

Pco Uco A 

A Phenolic nylon ablating, — — jr = 0.04 

Pco u co A 


Calculated 
Inviscid 


— — — — jr = 0.04, present method 
Pa} u co A 



(a) Hemisphere-cylinder, Re^ ^ = 9000. 


Measured 
O Nonablating 

□ Teflon ablating, — ® — r = O.OI 
Pco u co A 

O Delrin abloting, — — x- =0.08 

Pco u co A 

A Phenolic nylon ablating, — — -r = 0.04 

Pco u ao A 

Calculated 


Inviscid 



(c) Cone-cylinder, Re^ j = 9000. 


Measured 
O Nonablating 

□ Teflon ablating, m - A - = 0.24 

pco u co A 

O Delrin ablating, -— 77 — 7 r = 0.44 

PcO U G0 A 

A Phenolic nylon ablating A — tt = 0.18 

pco uco A 

Calculated 
Inviscid 


— 7 = 0.44, present method 

Pco u CO A 



(b) Hemisphere-cylinder, Re ro>c j = 1000. 


Measured 
O Nonablating 

□ Teflon ablating, m — r = 0.01 

Pco Uco A 

O Delrin ablating, m • ~r = 0.22 

Pco Uco A 

A Phenolic nylon ablating, — ~ — r = 0.1 1 
Pco uco A 

Calculated 
Inviscid 



(d) Cone- cylinder, Re WiC | = 1000. 


Figure 6.- The surface pressures for the bodies with ablating noses; M ro = 14. 
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The pressure data for the bodies at - 9000 (figs. 6(a) and 6(c)) 

show essentially no effect of nose ablation wliere measurements were made 
(beyond two diameters downstream of the nose) . The theory for the 
hemisphere-cylinder agrees with this . 

The data for the hemisphere-cylinder at Re^d = 1000 (fig. 6(b)) show 
increased pressures due to ablation just downstream of the nose, but the 
pressures rapidly approach the nonablating values farther downstream. The 
cone-cylinder data (fig. 6(d)) also are the same as the nonablating pressures 
far downstream of the nose, the only region where measurements with ablation 
were made. The viscous-interaction calculations for the hemisphere-cylinder 
agree qualitatively with the measured trend, but underpredict the level of 
the pressures for the Delrin nose. To see if this disagreement was due to 
blunting of the Delrin hemisphere, a test was made with a nonablating nose 
having the final shape of the ablated Delrin nose. The measured pressures 
agreed with the nonablating hemisphere pressures, so the disagreement cannot 
be attributed to nose blunting. 

It thus appears that nose ablation can locally increase surface 
pressures just downstream of the nose, but the effect does not persist very 
far downstream. This conclusion would not necessarily apply to entirely 
ab lating bodies . 


CONCLUSIONS 


From the results of the experiments and comparisons with viscous - 
interaction theories, the following conclusions were reached: 

1. The shock wave for the nonablating hemisphere-cylinder agreed with 
the predicted inviscid shock wave in both the boundary -layer and merged-layer 
regimes. The shock wave for the nonablating cone-cylinder was displaced out- 
ward from the predicted inviscid shock wave an increasing amount with 
decreasing Reynolds number over the range tested. Nose ablation displaced 
the shock waves outward for both bodies. All results were adequately 
predicted by viscous-interaction theory. 

2. The measured surface pressures on the cylinder were greater than the 
calculated inviscid pressures near the shoulder for both nonablating noses in 
the boundary- layer regime, but the differences rapidly diminished farther 
downstream. The results were adequately predicted by either of two viscous - 
interaction theories. In the merged-layer regime, measured pressures fell 
much below the predictions of the viscous-interaction theories. This result 
appears to be due to the merging of the shock wave and boundary layer within 
the shock layer causing a departure from a Rankine-Hugoniot shock wave. It 
thus appears that because of merging, highly cooled blunt bodies, at least 
for the first few nose diameters downstream, will never show the large 
viscous induced increases in surface pressure that have been measured for 
adiabatic-walled blunt bodies at higher Reynolds number. 


10 


3. Nose ablation in some cases increased cylinder pressures just 
downstream of the nose, but the effect does not persist very far downstream. 
These trends were qualitatively predicted by a viscous-interaction method. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, Nov. 26, 1969 
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